/***************************************************************************
 *   Copyright (C) 2007 by Alexandre Masselot,,,   *
 *   alex@alex-laptop   *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU Library General Public License as       *
 *   published by the Free Software Foundation; either version 2 of the    *
 *   License, or (at your option) any later version.                       *
 *                                                                         *
 *   This program is distributed in the hope that it will be useful,       *
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of        *
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         *
 *   GNU General Public License for more details.                          *
 *                                                                         *
 *   You should have received a copy of the GNU Library General Public     *
 *   License along with this program; if not, write to the                 *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/
#ifndef BOLDELTAMASSPPM_H
#define BOLDELTAMASSPPM_H

#include "deltamass.h"
#include <math.h>

#ifndef INFINITY  // This is a GNU C extension. MSCV lacks it.
#include <limits>
#define INFINITY  std::numeric_limits<double>::infinity() 
#endif

namespace bol {

/**
  ppm based mass error

  @author Alexandre Masselot,,, <alex@alex-laptop>
 */
  class DeltaMassPPM : public DeltaMass
  {
    private:
      double m_tolrel;
    public:
      DeltaMassPPM():DeltaMass(),m_tolrel(0){
      }

      ~DeltaMassPPM(){}


      bool isInRange(const double& m1, const double& m2)const{
        return (fabs(m1-m2)/(m1+m2))<=m_tolrel;
      }
      /**
       * the range of tol ppm around a mass
        *  warning, as we enforce a symmetric "distance" in ppm,
        *  sup=m*(1+m_tolrel)/(1-m_tolrel);
        *  inf=m*(1-m_tolrel)/(1+m_tolrel);

       */
      void range(const double& m, double& inf, double &sup)const{
        if(m_tolrel==0){
            inf=sup=m;
            return;
        }
        if(m_tolrel>1){
          sup=INFINITY;
          inf=0;
	  return;
        }
        sup=m*(1+m_tolrel)/(1-m_tolrel);
        inf=m*(1-m_tolrel)/(1+m_tolrel);
      }

      void setTol(double theValue) {
        m_tolrel = theValue*0.5e-6;
      }
  

      double tol() const {
        return m_tolrel*2.0e+6;
      }

      string toString()const{
        return boost::lexical_cast<string>(tol())+"ppm";
      }


  };

}

#endif
